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Aims. Mathematical modeling of bodily tides can be carried out in various ways. Most straightforward is the method of complex 
amplitudes, which is often used in the planetary science. Another method, employed both in planetary science and in astrophysics, 
is based on decomposition of each harmonic of the tide into two bulges oriented orthogonally to one another. We prove that the two 
methods are equivalent. Specifically, we demonstrate that the two-bulge method is not a separate approximation, but ensues directly 
from the Fourier expansion of a linear tidal theory equipped with an arbitrary rheological model involving a departure from elasticity. 
Methods. To this end, we use the most general mathematical formalism applicable to linear bodily tides. To express the tidal amend- 
ment to the potential of the perturbed primary, we act on the tide-raising potential of the perturbing secondary with a convolution 
operator. 

Results. This enables us to interconnect a complex Fourier component of the tidally generated potential of the perturbed primary with 
the appropriate complex Fourier component of the tide-raising potential of the secondary. Then we demonstrate how this interrelation 
entails the two-bulge description. 

Conclusions. While less economical mathematically, the two-bulge approach has a good illustrative power, and may be employed 
pH . on a par with a more concise method of complex amplitudes. At the same time, there exist situations where the two-bulge method 

p») 1 becomes more practical for technical calculations. 
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1. Introduction and aim 

On several occasions, it was suggested by different authors to model bodily tides with superposition of two symmetrical bulges. One 
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bulge is always aimed at the secondary, and thus implements the instantaneous reaction of the primary's shape and potential to the 
tide-rising gravitational pull exerted by the secondary. This portion of the tide is called "adiabatic tide" (Zahn 1966a,b) or "elastic 
tide" (Ferraz Mello 2012; Krasinsky 2006). The second bulge is assumed to align orthogonally to the direction to the tide-raising 
J^-j' | secondary, and thus is set to implement the entire nonelastic portion of the primary's deformation. This, second bulge is called 
"dissipative tide" (Zahn 1966a,b; Krasinsky 2006) or "creep tide" (Ferraz Mello 2012). 

In this note, we demonstrate that the two-bulge method is not a separate approximation, but ensues directly from the Fourier 



o> 



o 



, expansion of a linear tidal theory equipped with an arbitrary rheological model involving a departure from elasticity. 
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2. The static linear theory of bodily tides 

Let a spherical primary of radius R be subject to the gravitational pull by a secondary of mass M sec , residing at the position 
r * = (r*, (f>* ', A*) , where r* > R . At a surface point R = (R, (f>, A) of the primary body, the tidal potential generated by the secondary 
can be expanded over the Legendre polynomials P;(cos y) as 

W(R,r*) = V W,(R, O = V ( — ) P,(cosy) . (1) 

Here G denotes Newton's gravity constant, while y is the angular separation between the vectors r * and R pointing from the 
primary's centre. The longitudes A, A* are measured from a fixed meridian on the primary body, the latitudes 0, (p* being reckoned 
from the equator. The index / is conventionally named as the degree. In (Q]i the / = term is missing, because it corresponds to the 
principal, Newtonian part of the secondary's potential, and is not a part of the perturbation. Omission of the / = 1 term is a more 
subtle point related to the fact that we are describing the motion of the secondary relative to the primary body, and not relative to an 
inertial frame (see equations 8 - 1 1 in Efroimsky & Williams 2009). 

Within the linear theory, the / th term Wi(R , r *) of the secondary's potential generates a linear alteration of the primary's shape. 
This alteration, in its turn, causes a linear amendment Ui(r) to the gravitational potential of the primary, where linear means: linear 
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in W/(R , r *) . The theory of potential requires that outside the primary body Ui(r) should scale with the distance as r </+;) . Hence 
the said change in the primary's potential may be written down as 

Ui(r) =**(-) Wi(R, r ') , (2) 

where R is the mean equatorial radius of the primary, R = (R , <p , A) is a point on the primary's surface, while r - (r , <p , A) 
is an exterior point right above the surface point R , at a radius r > R. The numerical factors k\ are the degree-Z Love numbers 
calculated from the rheology of the primary body. 

The overall tidally caused change of the primary's potential thus amounts to 

U(r) = J] Ufa) = J] k, - WAR, r*) . (3) 

1=2 1=2 r 

3. Dynamical linear theories of bodily tides 

In realistic situations, the disturbing potential is evolving, so equation (Q]) assumes the form of 

i 

(4) 



W(R,r*(t)) = V W,(R, r*(t)) = V — P,(cos r (f)) 



/=2 v 7 1=2 

Then one should expect the distortion of the primary, as well as the corresponding amendment to its potential at an exterior point 
r , to become a function of time: U (r, f) . 

3.1. Elastic dynamical tides 

Had the tides contained only instantaneous, elastic components, the expressions for the tidal potential would mimic (j2]-|3]l. At each 
instant of time t , the degree-Z term of the tide-raising potential W of the orbiting secondary would generate instantaneously an 
appropriate degree-Z term of the tidal potential of the primary: 

(elastic) I R 

Ufa, t) = kA - J WAR, r\t)) , (5) 
so the total tidal amendment to the potential of the primary would look: 

DO DO n 7j-1 

(elastic) (elastic) ^ / R\ 

U(r, = 2j Ul{r > f) = Z kl (-) Wl(R > r {t)) ■ (6) 

1=2 1=2 r 

Needless to say, in realistic materials the internal friction prevents the tidal deformation from being instantaneous. 

3.2. Realistic dynamical tides 

To describe deviation from elasticity, we spell the two basic assumptions whereon a linear dynamical theory of bodily tides is based: 
[1] tidal deformation is linear with respect to the stress generated by the tide-raising potential; 

[2] the deformation is not fully elastic: it incorporates both an immediate and delayed portions (delayed - relative to the tide-raising 
potential). 

Mathematically, assumption [1] means that an infinitesimal increment AW/(ff , r *(/')) , whereby the perturbing potential in- 
creased at the time t ' in the past, results in a proportional present-time increment of the tidally distorted shape of the primary and, 
accordingly, in a proportional increment of the tidal amendment to its potential: 



MJAr, f) = (-) *,(*-*') AWAR, r *(/')) , (7) 

kfa - t ') being a function describing the delayed reaction of the shape. The mathematical form of this function is defined by the 
rheology of the body and by its self-gravitation. 

Thus the addition Ui to the primary's potential gets expressed through the tide-raising potential Wi by a linear integral operator: 

UAr, = ^ kAt -t')Wi{R,r \t '))dt' , (8) 



overdot denoting a time derivative. Integration of d8) by parts renders: 
Ufa, t) = [kAOWfa - kA°°WA-°°)] + j kAt-t') WAR,r*(t'))dt' , 



l+l 

(9) 
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where the relaxed term - A:/(oo)W(-oo) should be neglected. Indeed, the current events cannot be influenced by the perturbation 
Wi(-ca) in the infinite past, wherefore A:/(oo) = . Of the remaining two terms, the unrelaxed term ki{Q)Wi(t) reflects the elastic 
part of the deformation, while the integral expresses the delayed components - viscous and anelastic. 

Assumption [2] means that both the unrelaxed term ki{Q)W/(t) and the delayed term given by the integral should be kept. As 
explained in Efroimsky (2012a,b), the relaxed term may be easily incorporated into the integral, where it should show up multiplied 
with a Heaviside step function ®(f - t ') . Then our expression for the tidal potential will acquire the simple form of 

U,(r, t) = J k,(t-t')W,(R, r*(t'))dt' , (10) 

k[(t-t ') now incorporating both the delayed-reaction terms and the elastic term ki(0) 0(? — t'). The elastic part will enter the kernel 
kAt—t ') as ki{0) 6{t -t'), with 5{t — t') being the Dirac delta function. Integration of this term will furnish fc/(0)W/(f) , as in ©■ 



4. Fourier components of tidal stresses and strains 

4. 1 . Tidal modes and tidal frequencies 

The sidereal angle and the spin rate of a tidally-perturbed primary are normally denoted with 8 and 8 , while the node, pericentre, 
and mean anomaly of a tide-raising secondary, as seen from the primary, are denoted with O , u> , and At. 

In the Darwin-Kaula theory, the tide-raising potential W, the primary's deformation, and the tidal amendment U to the primary's 
potential are expanded over the modes 

oJimpq = (I ~ 2p) ti> + (I -2p + q) M + m (O - 6) « (/ - 2p + q) n - m 8 , (11) 

/, m, p, q being integers, and n being the mean motion. Dependent upon the values of the mean motion, spin rate, and the indices, 
the tidal modes u>i mpq may be positive or negative or zero. 

The actual forcing frequencies of the resulting stresses and strains in the primary's material are the absolute values of the tidal 
modes: 

Xlmpq = I Wlmpq I , (12) 

so these frequencies are always positive. 



4.2. Fourier expansions 

In practical calculations, it is extremely convenient to employ complex stresses and strains, under the convention that the actual, 
physical quantities are the real parts of their complex counterparts. This way, the Fourier series for the stress and strain look: 

oo oo 

o-yvif) = ^ o-yy&s) cos [x s t + <Po-(x s ) ] = J] 1<e [ o-yviXs) exp(i^)] , (13) 

CO CO 

Uyvit) = ^ Uyy(x s ) co$,\x s t + tp u (x s )\ = ^^e[s yv O s ) exp(ix s t) ] , (14) 
yv being tensor indices, and s being a concise notation for Impq. The complex amplitudes are 

o-yvix) = o- yv (x) exp [iVo-fr)] , u yv (x) = u yv (x) exp[i(p u (x)] , (15) 

where the initial phases <ficr(x) an d <Pu(X) are chosen so that the real amplitudes o- yv {x s ) and u yv (x s ) are non-negative. 
For a continuous spectrum, the sums get replaced with integrals over frequency: 



(16) 



'X' \ _ 'X' 

Cyvix) e dx and u yv (t) = %e I u yv (x) e dx ■ 
o Jo 

Similarly, the tide-raising potential Wi and the potential Ui of the primary get expanded into a sum or integral over the tidal modes: 

XOO . , /-»O0 . , 

W,(u})e do> and U,(t) = He U,(co) e du , (17) 
CO xj — OO 

where the complex amplitudes are expressed via the real amplitudes and the initial phases by 

Wi(oj) = wm exp[i<p Wi ((o)] , Um = U,(oj) expfiV^M] . (18) 

The phases <p w (co) and 'P Ui {cl>) can always be set in such a way that the real amplitudes Wi(x) and Ui(x) are non-negative. 

Both in ( fTol and ( fTTl i. the actual, physical spectral components are the real parts of the complex ones. However, there also is an 
important difference between (fTol l and ( fTTl i. While the stresses and strains are habitually expanded in ( fTol l over positive frequencies 
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X only, the potentials in ( fTTl i are expanded over the tidal modes to, which can be positive or negative or zero, as demonstrated in 
the Darwin-Kaula theory of tides. 

It is of course a common fact that a real function can be decomposed into a Fourier series or integral over only positive frequen- 
cies. This way, expansions of the potentials over x — M ma y appear to be sufficient, because a contribution from some negative 
tidal mode to < can be shown to coincide with the contribution from the appropriate positive mode \u>\ > . Therefore, (TPTl i may 
be rewritten simply as 

Wi(x)e' X d X and U t (t) = <Re Ufa)e' X d X , (19) 
o Jo 

where Wfa) = 2 Wi(cj) and Ufa) = 2 Ufa) . 

Surprisingly, the theory of tides is a rare exception from the rule, in that this theory does distinguish between the contribution 
from a negative tidal mode and that from a positive mode of the same absolute value. Fortunately, the difference shows up only at 
the stage when one calculates tidal forces or torques (Efroimsky 2012a,b). As in the current paper we discuss potentials only, we 
shall ignore this subtlety and shall employ ( fT9l instead of (fTTl i. 

4.3. Dynamical analogues to the Love number 
Insertion of ( fT9l into ([Tol l entails 

Ufa) = 1 fa) Wfa) (20) 

or, in a more detailed form: 

\Ufa)\e = (-) \kfa)\ e \Wfa)\ e , (21) 

where the complex function 

kfa) = \kfa)\e~' eM , (22) 

is a Fourier component of the kernel k t {t - t ') of the integral operator (fLOt . The kernels k/(t — f ') are named Love functions, a 
term suggested by Churkin (1998). The functions kfa) are called complex Love numbers, in understanding however that these 
"numbers" change with frequency. An approach similar to (l22l i was taken by Mathis and Le Poncin Lafitte (2009). These authors 
introduced a complex impedance as the ratio between the complex Love number and the static Love number. Their equation (107) 
is equivalent to our (l22l . up to a convention on the sign of the argument. 

We see from (ETT i that at each frequency x > the negative argument efa) is a measure of lagging of the spectral component 
Ufa) , relative to appropriate spectral component Wfa) : 

Vu.Od = <fiw,M - efa) . (23) 

5. Decomposition of a tidal mode into an in-phase part and an in-quadrature (lagging by 90°) part 

In the expression d2TT >. we may set the initial phase of Wfa) to be zero, and reckon the phase of Ufa) from Wfa) . This will enable 
us to single out, in Ufa) , a part which is in phase with Wfa) , and also to see what part of Ufa) is out of phase with Wfa) . 

According to (1231 1. nullification of the phase of Wfa) makes the phase of Ufa) equal to the negative argument of kfa) . So 
(f2Tb will assume the form of 

U h fa)e = (-) k lQ (x) e W h fa) e , (24) 

where we introduced simplified notations Ui (x) , Wi (x) and ki for real amplitudes: 

Ui Q (x) = Wfa)\ , W lQ (x) = \ Wfa)\ , k l0 (x) = \kfa)\ ■ (25) 
By means of the Euler formula, (l24l can be trivially expanded as 
Ui (x) [ cos (xt - efa) ) + i sin (xt - efa) ) ] 



2/+1 



= (-) ^ W [cos efa) -i sine/Or)] W/ tr) [ cos(#f) + i sinOff) ] . (26) 

The actual, physical tide-raising potential Wfa) is the real part of the complex Wfa). So it is rendered by Wnfa) cos(^f) , where 
we set the initial phase nil, as agreed above. Similarly, the actual, physical tidal potential Ufa) is the real part of the complex Ufa), 
and it reads as: 

Ui (x) cos efa) cos(xt)- U K fa) sin efa) sin(xt) 

= (-) k u fa) cos efa) W, Q (x) cosOrf) + (-) k K fa) sinefa) W, ( fa) sin(xt) . (27) 
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In this expression for the actual real U\ , we see not one but two terms. One is the elastic part of the tide, a part that is in phase with 
the real Wi . This is the term proportional to cosfj^f) : 

(in phase) {R\ 2l+l 

Ui = Uio(x) cos ei(x) cosCff) = ^-J k IQ (x) cos q(x) W, (x) cos(yf) , (28) 

We see that the "response factor" is equal to k^ix) cos e/(x) , where ki (x) is the real amplitude of the complex Love number (call 
this amplitude dynamical Love number). 

The second component is the in-quadrature term, which is proportional to sin(yf) : 

(in quadrature) I R\ 2l+l 

Ui = - U l0 (x) sine/Or) sin^f) = [-) k, Q (x) sin e t (x) Wi (x) sinCyf) . (29) 

Naturally, the expression of this term via Wi contains ki (x) sin e/(x) . 

Therefore, as soon as we express Ui via Wi by an integral operator permitting delayed action (Eqn. 9), we automatically arrive 
at the two components of the tidal {//, at each frequency involved, - the adiabatic component and the dissipative component (using 
the terms of Zahn 1966) or the elastic one and the creep one (as Ferraz-Mello 2012 named them)Q 

Similar decomposition will take place for the harmonic modes of the primary's surface elevation, except that the Love number 
hi will be involved instead of k/ . For the principal, {Impq} = {2200} , tidal mode, this situation is illustrated, in a very exaggerated 
manner, by Figure 1 . 




e 

Fig. 1. Decomposition of the semidiurnal tide into the elastic and dissipative components. 

The primary is spinning at the rate 8 , with 8 being its sidereal angle. The secondary is orbiting the primary with the mean motion n . 
The ellipse shaded green depicts the principal, semidiurnal mode of the tides exerted on the primary by the secondary. This is the mode 
numbered with {Impq} = {2200). The semimajor axis of the green ellipse deviates from the direction towards the perturber by a ge- 
ometric lag angle. For the semidiurnal tide, this angle is equal to j 62(^2200) > where 62 is the phase lag as a function of the prin- 
cipal tidal mode W2200 • 

It is shown here, in an exaggerated manner, how the principal tidal mode gets decomposed into two components. One is the elastic bulge de- 
picted by the ellipse aligned with the perturber. Another is the dissipative bulge depicted with the ellipse aligned orthogonally. 



It should finally be mentioned that in the case of stars and giant planets it is, technically, difficult to solve the emerging hydro- 
dynamical problem analytically. Therefore, in practical calculations of fluid equilibrium tides in such bodies, the two-bulge method 
is the only known elegant option of getting an acceptable analytical solution. One first solves the hydrodynamical equations gov- 
erning the adiabatic component d26l) . The adiabatic adjustment of the structure is deduced from the hydrostatic balance. Then, since 
the velocity field is preserved along isobars, one obtains the radial component of the adiabatic velocity field induced by the tide. 
The subsequent calculation of the horizontal component is based on the velocity field being divergence-free. The knowledge of the 
adiabatic velocity field makes it possible to determine the viscous force. This force drives the dissipative component of the velocity 
field. The process leads to mass redistribution inside the body, and thus to perturbation of the density and gravitational potential 
d27b . The described, two-bulge approach is available because in the fluid equilibrium tide the dissipative component is much weaker 
than the adiabatic component. The method is implemented, e.g., in the work by Remus et al. (2012a), which furthers the original 
two-bulge approach offered by Zahn (1966). 

1 Also see the work by Krasinsky (2006), who employed the terms elastic and dissipative. 
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6. Conclusions 

In this short note, we pointed out a simple rule linking the two methods (or, possibly better to say, two languages), in which bodily 
tides have been described by different authors. The language of complex amplitudes is more economical and is conventional to 
those who studied the theory of vibrations, in physics or engineering. The language of two bulges turns out to be equivalent to that 
of complex amplitudes. Although less economical mathematically, the two-bulges language has some illustrative power, and may 
be employed (like, e.g., in Remus et al. 2012a,b) on a par with the more concise method of complex amplitudes. 

Importantly, the existence of two mutually orthogonal bulges at each tidal frequency is not a separate approximation (as was 
presumed by some of the devotees of the two-bulge method), but is a consequence of the linearity assumption implemented by 
the integral operator (ITOV As soon as we say that the tidal deformation is linear but not fully elastic - this deformation can be 
decomposed, at each tidal frequency, into an in-phase and an in-quadrature part. 
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